Equilibration of a hard-disks system 



V.M. Somsikov 

Laboratory of Physics of the geogeliocosmic relation, 
Institute of Ionosphere, Almaty, 480020, Kazakhstan 

Abstract 

The process of equilibration of a colliding hard-disks system is studied in the 
framework of classical mechanic. The method consists of dividing the nonequi- 
librium system into the interacting subsystems; the evolution one of these sub- 
systems is analyzed employing generalized Lagrange and Liouville equations. The 
subsystem-subsystem interaction force is considered as an evolution parameter. 
The mechanism by which its system equilibrates is described. 

Keywords: nonequilibrium; irreversibility; many-body systems; entropy; evolu- 
tion. 



1. Introduction 

The existing methods of analyzing many-body systems based on the models of disks and 
spheres have been shown to be highly effective, see e.g. [Kryilov, 1950; Sinai, 1970; Beijeren & 
Dorfman, 1995]. Nevertheless, it is not clear whether it is possible to prove rigorously, within 
the framework of classical mechanics laws, that equilibrium is established in a sum system, 
see [Petrosky &: Prigogine, 1997; Zaslavsky, 1984,1999]. Such an enterprise would demand new 
methods and new approaches. 

Some features of the nonequilibrium evolution of colliding disks were studied and reported in 
[Somsikov, 1996, 1998, 2001a]. This approach enables one to carry out analytical and numerical 
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study of the evolution process within the framework of classical mechanics on the basis of 
equation of motion of hard disks. These equations are written down for a matrix of colliding 
disks. It was found that the establishment of equilibrium was determined by the rate of decrease 
of forces between disks group, [Somsikov, 1996,2001b]. 

In the present work the process of establishment of equilibrium state for a hard-disks system 
within the framework of the laws of the classical mechanics is studied. The principal feature of 
this study is that the interaction force of subsystems into which a disks-system is divided, is 
used as evolution parameter. 

This allows us to avoid application of probabilistic concepts for an explanation of the mech- 
anism of equilibration. Moreover it helps us understand the nature of the probabilistic laws in 
a systems of the classical mechanics. 

The work is constructed as follows. Using D'Alambert principle and equations of motion for 
disks, we construct generalized Lagrange, Hamiltonian and Liouville equations for a subsystem, 
selected from the full system. With the help of these equations, we analyze peculiarities of 
evolution of the disks-system from a nonequilibrium state to equilibrium. 

The study is based on the following method. First we prepare a nonequilibrium disks system. 
A macroscopic subsystem of disks is selected from this system. The evolution of forces on this 
subsystem by other subsystems is analyzed. 

2. Generalized equations of the classical mechanics 

It is known that the methods, based on canonical Lagrange or Hamilton equations, can 
be used successfully for investigating many-body systems, which is not faraway from an equi- 
librium state see [Lanczos, 1962; Landau, 1976; Landau & Lifshits, 1974]. At the same time 
these methods have inherent problem when we attempt to use them for the study of strongly 
nonequilibrium systems, especially ones with nonholonomic connections and polygenic forces 
(see [Lanczos, 1962]). But in the real world, the majority of systems are of this type. Indeed, 
even a disks-systems in the space with ideally reflecting walls are nonholonomic [Goldstein, 
1975]. Therefore for considering the process of evolution to equilibrium state in a hard-disks 
system, we should modify the equations of the classical mechanic so that they are applicable 
to nonequilibrium systems. For deriving these equations, we must take into account the fact 
that the forces acting between subsystems, in general, cannot be expressed through a potential 
function; the work done in going from one point to the other in a configuration space depends 
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on the path taken. These equations wih allow us to consider process equilibration, by using the 
interaction force of subsystems, as a dynamic evolutionary parameter. In accord with Lanczos 
[1962], we shall name these as "generalized equations". 

Let us take a system, which consist of N disks. Divide this system into R subsystems, so 
that in each subsystem will be of T disks. Therefore, N = RT. The energy of the system is 
constant. It is equal to the sum of internal energies of all subsystems and interaction energies 
between subsystems. Let us selected one of them, which we call p-subsystem. Let 6W^, be the 
virtual work of active forces in p-subsystem. In the general case, this work can be express as 

T N-T T 

follows: dW^ = J2 J2 ^ks^fk = ^k^fk^ where k = 1,2, 3.. .T are disks of p - subsystem, 

k=l s=l k=l 

s = 1,2,3. ..N — T are external disks that interacts with the p -disks. F^^ is the interaction force 

N-T 

of k - and s - disks, Srk is the virtual displacement of the k -disk, = J2 ^ks- '^^^ virtual 

s=l 

work of the interaction force of internal disks of p-subsystem is equal to zero. 

In case of the pair interaction the virtual work of external forces is of the form: 5W^ = 

T T 

F^Srk = J2 ^ks^'^'k-, because in this case F'^ = Fj^^. The inertial force can be presented as 

k=l k=l 
T . 

^^in ~ S ^k^fk- The sum of the active and the inertial forces constitute as effective forces. 
k=l 

The principle of D'Alambert asserts that the work of effective forces is always zero [Lanczos, 
1962], i.e. 

6Wl = SW[^ - 8Wl (1) 

The feature of the virtual work above implies, that in general it does not get reduced to complete 
differential. 

R T 

Prom the equations of motion for hard disks see [Somsikov, 2001b], it follows: Yl, ( ^ks)^''^k = 

p=l k=l 

0. Therefore total active and inertial work for all subsystems at any moment of time is equal 

R R 

to zero, i.e. J2 ^^fn — S <^^a = 0. This equality can take place in two cases: when the sum 

p=i p=i 

of nonzero terms is equal to zero, or when each term of the sum is equal to zero. Obviously the 
second case, appropriate for an equilibrium state, takes place when T ^ oo. For this case, with 
the help of the equations of motion for hard disks [Somsikov, 2001b], it is possible to record: 

T T T 

E ^ks5{^Pks{t))/^ks{t)\/^ks{t)\ = J2 PL = E = 0' (2) 

k=l k=l k=l 

where Vk is a /c-disks velocity; cpks = i/Se^^'"'; i is an imaginary unit; j3 = siwdhs'- i-'^ks is 

scattering angle for k and s colliding disks, which varies from to tt; 6{x) is a delta function; 

t 

ipks = 1 — l^fesl; hsit) = ^ksdt are distances between centers of colliding disks; r^^ = 'r'k~f^ 
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are initial values of disks coordinates; A^g = ~ sue relative velocities. The impact 
parameter dks = cos'&ks is determined by the distance between centers of the colliding disks in 
the complex plane with real axis x and imaginary axis y. The value of the impact parameter 
is: dhs = Ii^ihs^ks) /\hs^ks\- The A;-disk swoops on the j- disk lying along the x - axis. Mass 
and diameter of each disk is set equal to 1. Boundary conditions are either periodical or hard 
walls. Equating the right-hand side of Eq. (2) to zero implies, that the selected p-subsystem is 
in a stationary state. 

To derive the general Lagrange equation for p-subsystem, let us transform D'Alambert 
equation (1) by multiplying it by dt, and integrating it over an interval from ti to t-2. In general 
we have, 

/ = / E [iVk -TFl-- Fl]5rudt = 

ti ti k=l j^k 

sJki ^kdt - ni{Fl + E F!k)Srk]dt - [E VkSvkWll (2a) 
ti k=l ti k=l jj^k fe=l 

T 

In Eq. (2a) the term E FL determines the force of interaction in the p-subsystem, k and 

T 

j are the colliding disks of the p-subsystem. The term E F^ is the force on the p-subsystem 

k=l 

from the rest. On demanded at the ends of the interval [ti, t2], the virtual displacements is zero. 
Then the last term in (2a) will be zero. 

Let us assume that a subsystem has equilibrated. Then for internal forces of interaction 

in a subsystem we can set a function U{ri,r2, ■■■Tt), for which the following condition is sat- 

t'z T T t2 

isfied: /[E J2 Fj^j6rii.]dt = —5 J U{ri,r2, ...rT)dt. Here ri,r2---rT are the coordinates of the 

tl k = l jy^k tl 

p-subsystem disks. In the general case it is impossible to express forces on p-subsystem, as a 
gradient see [Lanczos, 1962]. In this case the Eq. (2a), can be written as 

t2 t2 rp 



In Eq. (3) we denote i^p = E + U{ri,r2, .-rT)- If the interaction of disks is potential, then 

fc=i 

Lp will include also internal potential energy of the p-subsystem - J7(ri, r2, ...rr). Since any 
variations of integral in equation (3) will be zero,we can set: 

T 

In the above T, F^ = Fp. 
k=l 
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Eq. (4) is a generalized equation of Lagrange for a p-subsystem. Fp is the polygenic force 
acting on the p-subsystem. When Fp = 0, Eq. (4) transforms to a canonical equation of 
Lagrange for equilibrium, conservative system. The equality, Fp = 0, is a sufficient condition 
for a stationary state. 



Let us now derive Hamilton equation for the chosen p-subsystem. The differential for Lp can 



T 

be written as, dLp = J2 ("gf^t^^'fc + ov'dVk) + -gf-dt, where = pk is disks momentum. With 



fe=i 

T dL 

the help of Lagrange transformation, it is possible to write: d[ J2 Pk^k — Lp] = (— -g^f^^fe + 

k=i k=i 

T 

Vkdpk) Qt'd^- Since -gf- = — q^, where Hp = J2 Pk^k ~ -^pj we have from (4), 

fe=i 



dr. 



= -A + ft- (5) 

fi = Vi. (6) 

dpk 

The above is a general Hamilton equations for the selected p-subsystem. The right-hand side 
of Eq. (5) denote the external forces, which act on p-subsystem. 

Using Eqs. (5) and (6), we can find the Liouville equation for p-subsystem. For this 
purpose, let us take a generalized current vector - Jp = {r'k,Pk) of the p-subsystem in a phase 
space [Zaslavsky, 1984]. Prom Eqs. (5) and (6), we find: 

The differential form of particles number conservation law in the subsystem is a continuity 
equation: ^ -|- div{Jpfp) = 0, where fp = fp{rk,Pk,t) is a normalized distribution function of 
disks in the p-subsystem. With the help of continuity equation and Eq. (7) for divergence of a 

T 

generalized current vector in a phase space, we can show that: ^ = ^ -|- Yl (^k^+Pk§^) = 

T 

^ -I- div{Jpfp) - fpdivJp = - fp Y Jr^fe • Thus, we have: 

fe=i 

Equation (8) is a Liouville equation for p-subsystem. It has a formal solution: 

fp = const • exp [- / ( E J^^D^*] (8a) 
fe=i 

Prom this solution it follows, that the p-subsystem will be in a stationary state when the external 

T 

forces disappear, i.e./ ( Y ^^Fl^)dt = 0. So, the exponent index of the solution of the Eq. (8) 

k=l 



determines the characteristic relaxation time of the system to the equiUbrium state. Hence, 
the change of the phase volume of the p-subsystem will last during time of aspiration of the 
force, Fp, to zero. How this comes about is discussed bellow. It is possible for all points of the 
phase space except for "islands", filled by periodic and quasiperiodic orbits of some Hamiltonian 
systems, see [Loskytov & Mihailov, 1990; Zaslavsky, 1999]. If at the time of preparation, the 
system is in such an island, it will return to it periodically. For periodic points wc do not have 
mixing, and the correlations do not disappear. We can then say that the probability of return 
to the initial periodic point is determined by the probability of preparing the system at its 
periodic point. 

Let us consider the question: how is the description of selected subsystems connected to 
disks system description as a whole. 

R T 

As, Yl S -^fc = 0) the next equation for the full system Lagrangian, Lr, will be: 

P=i fe=i 

dt dVk dru ^ ' 

and the appropriate Liouville equation is 

dfR . dfn . djR 
dt dvk opk 

The function, /r, corresponds to the full system that is conservative. Therefore, wc have: 

R J J ^ 

divJp = 0. This expression is equivalent to the next equality: ^( In/p) = ^(1^ H fp) = 
p=i p=i p=i 

R ^ R R R 

in fp) ^( n /p) = 0- So, Yl fp = const. In equilibrium state we have Yl fp = fR- Because 
p=\ p=\ p=\ p=\ 

R R 
the equality Fp = is fulfilled at all times, we have the equality, 11/3? = Z^' ^ integral 

p=i p=i 
of motion. This is in agreement with Liouville theorem eabout conservation of phase space 

[Landau & Lifshits, 1973]. 

The simultaneous fulfillment of conditions of phase volume preservation for full system and 

validation of time-dependence solution (b) of Eq. (8) for nonequilibrium subsystems is correct 

only when the centre of mass of a system moves on trajectory, reversible in time. It should 

take place irrespective of, whether all subsystems are in equilibrium or not. Therefore, for the 

considered p-subsystem only those irreversible redistribution of phase volume and energy that 

lead to reversibility of motion of the centre of mass of all system, is possible. 

3. The mechanism of equilibration 
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Let us show that for a hard-disks system, the external force decreases due to the mixing 

T 

properties, i.e. / ( l^^k)^^ ~^ ^ when t ^ oo. This means that because of mixing, the 

k=i ^'^ 

system will go to a stationary state irrespective of where it was at the initial time (except the 
periodic points). 

The mixing properties of two hard disks was proved by Sinai [1970] and Zaslavsky [1984]. It 
was shown that the correlation function for colliding disks is: R{t) exp{nlnK), where K = p/2: 
p-length of free run; n-the number of collisions. Therefore the characteristics time of the decay 
of correlations, td, for regular collision disks with unit diameter through a time interval r is 
determined by equation: td = t/ {In p/2). The condition K = p/2 > 1 is satisfied by rare 
gas. Hence, two hard disks are mixed. As will be shown bellow, the mixing property is both, 
necessary and sufficient for the equilibration. Although in the situation when we have several 
disks, the strict mathematical proof of a mixing property is absent; nevertheless the existence 
of the mixing property is usually accepted a priori as well. 

For proving the property of aspiration of the resulting force to zero for a hard-disks system, 
we shall simplify Eq. (2) for the p-subsystem. Let us assume, that all disks collide simultane- 
ously in equal, short enough intervals of time, r. It is clear, that if for this condition the system 
goes to equilibrium then does so definitely for the general case. After such simplification, Eq. 
(2) for a p-subsystem, can be written as 

= V\sin)^til) (11) 

Here, to each "A;" disk from jo-subsystem in moment of time, nr corresponds to "s(n)", disk 
from other subsystems; n = 1,2, 3, .... 

The evolution of the p-subsystem is determined by the vector, V^, with components denoting 
the velocities of disks of the p-subsystem: = {V^},k = 1,2,3, ...,r. Some of the time- 
dependent properties of this subsystem will be determined by studying the sum of it components. 
Let us designate this sum as Tp. Carrying out the summation in (11) over all disks of the p- 
subsystem, we obtain: 

t," = E^Vn)AL-(;.)=i^P (12) 
k=l 

Equation (12) describes the change of total momentum, acting on the p-subsystem as a 
result of collisions at time nr. The aspiration of a total momentum to zero is equivalent to 
aspiration to zero of the force, Fp. 
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Now let us show, if the mixing property for a disks system is assumed, the homogeneous 
distribution of impact parameters of disks occur as well. 

In accordance with the mixing condition, we have the following [Loskytov k, Mihailov, 1990], 

^i{6)/^i{d) = 6/d (13) 

Here, fi{d), is a measure corresponding to the total value of impact parameter - "d"; S is an 
arbitrary interval of the impact parameter and, fJ,{6), is a corresponding measure. Equation 
(13) implies the number of collisions is proportion to the interval "(5". It also implies that the 
distribution of the impact parameters is homogeneous. 

As is well known, see e.g. [Loskytov &: Mihailov, 1990; Zaslavsky, 1984], for mixed systems 
of correlations decay. For Eq. (12) this condition can be written down as < fks{n)^k^(n-i) 
'^ks{n) -^^^ ^fcs(n-i) ^ average from two multiplied functions is equal to the multiplica- 

tion of the average of these functions, (p'^g^^^ come from impact parameters, and come 
from relative velocities of colliding disks. Therefore this condition is similar to the condition of 
independence of coordinates and momenta, widely used in statistical physics see e.g. [Rumer 
& Ryvkin, 1977]. 

Thus, it is possible to carry out the summation in the multiplier, V'fes(n)' ^^^^ impact pa- 
rameters, independent of the summation of A^~^_^-^ over relative velocities of colliding disks. 
Then, under the condition of the homogeneous distribution of impact parameters and when 
r >> 0, we can go from summation to integration. We will have, see [Somsikov, 2001a]: 

rp TT 

= l/T^lim Y.^ks = ^j ^Ld{cos^) = -| (14) 
k=i 

where G = 2 is the normalization factor. 

Taking into account (14), we have from Eq. (12): 

2 ^ 

k=l 

The negative sign in the right-hand side Eq. (15) means, that the force, Fp, decreases. 

The next question is about the stability of a stationary point. The stability of a stationary 
point of p-subsystem can be established with the help of the Eq. (15). We set the initial 
deviation from stationary point, and then consider, how this deviation changes with time. 

Let the point, Zq, be a stationary point, so that, Fp, acting on p-subsystem, is zero. Prom 
the Lyapunov's theorem about stability it follows that the point, Zq, is asymptotically stable if 
any deviation from it gets attenuated. 
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Let us expand the left and right-hand sides of the equation (15) in a series by small pa- 
rameter, f, of perturbation of velocities of disks of the p-subsystem, near point, Zq, and keep 
terms up to first-order infinitesimal. The expansion of the left-hand side of the Eq. (15) gives: 

T 

V = £jk, where the summation is carried out on the components of the variation, v. In the 
k=l 

T 

expansion of the right-hand side, there remains only (— o) J2 = — 2/3u. A contribution into 

fe=i 

the expansion is given by collisions of disks of the p-subsystem, with disks of its complement. 
We have: 

= (16) 

Equation (16) means, that any deviation from an equilibrium state will decay. Hence, the 
stationary point in the presence of mixing is steady. Stability is provided by emergence of 
returning force, Fp, at a deviation of a subsystem from an equilibrium point. We shall note, 
that the Eq. (16) also follows from the theory of fluctuations see [Landau, 1976]. 

Let us show, that the emergency of force due to deviation of system from its equilibrium, 
provides restriction of spontaneous fluctuations. 

Take the system in a nonequilibrium condition. As follows from the previous statement, 
any nonequilibrium condition is characterized by the force, Fp, which acts on the p-subsystcm. 
This force is determined by Eq. (12). The time decrease of the force, Fp, is determined by 
equation: tdin = / Thus, if the system somehow appears at a nonequilibrium point, over a 
characteristic time, tdin, then it should return to equilibrium. 

For the further consideration we shall accept two statements which follow from the mecha- 
nism of equilibration. 

First: the degree nonequilibrium is defined by the force, Fp. So, there is a mutual unique 
conformity between Fp and a phases space points. 

Second: we shall consider, that the spontaneous deviations increase of system from an 
equilibrium condition occurs under increase of Fp. The deviation is proportional to force. 

If these conditions are fulfilled, it is possible to prove that such fiuctuations are realized only, 
if tdin > iprob- Time, tprob is determined by probabilistic principles. According to the formula 
Smoluhovsky [Zaslavsky, 1984], for the case of an ergodic system, average resetting time, t, 
or Poincare' cycle time, is equal to tprob = t{l — Po)/{Po — Pi), were Pi is the probability of 
reversibility during the time t, Pq is the probability of initial phase region. Suppose that the 
system in a probabilistic way begins to deviate from equilibrium. The characteristic time of a 
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deviation to any point, Zp, should be: ~ iprob- But because in this point the force, Fp, acts 
on subsystem as returning, the system will approach to the point, Zp, in a probabilistic way if 

only tprob < tdin- 

Thus, within the scope of assumption made in this work, the dynamics of a hard-disks sys- 
tem is completely determined by the deterministic classical mechanics. Therefore, the need for 

probabilistic principles in the description of evolution of the system and its area of their use are 
determined by roughness of transition from summation to integration on impact parameters, 
and also from the periodic or quasiperiodic points at the time of preparation. 

4. Conclusion 

Splitting conservative, nonequilibrium system as a set of interacting subsystems, and the 
analysis of the evolution of forces of interaction between these subsystems is a basic idea pro- 
posed in this paper. This approved allows us to take into account the exchange of energy 
between subsystems and a feedback between active and inertial forces of their interaction. The 
feedback provides relaxation of the system to equilibrium. Therefore this approach is applicable 
for studying evolutionary processes in open systems with polygenic forces and nonholonomic 
connections. 

The analysis nonequilibrium evolution is based on generalized equation of Lagrange set for 
a selected subsystem. The right-hand side of this equation is a polygenic force of interaction 
of subsystems. Because of mixing, this force tends to zero. It causes evolution of system to 
equilibrium state. The establishment of equilibrium is possible starting from all the points of 
the phase space except for periodic points. In the vicinity of equilibrium, where, Fp is near 
zero, the statistical theory of fluctuations constructed on the basis of the canonical Hamilton 
equations is applicable. 

Stability of an equilibrium state is ensured by the emergence of the returning force when a 
system deviates from equilibrium. It imposes appropriate restrictions on amplitudes of probable 
fluctuations of system. 

Within the framework of this study, the dynamics of the system is deterministic. Prob- 
abilistic principles enter only by the uncertainty of initial conditions, and coarse graining of 
transition from the discrete to the continuous. 

Let us compare our explanation of irreversibilities with existing one, see, for example [Za- 
slavsky, 1984]. In accordance with this explanation, irreversibility is a consequence of "coarse 
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graining" of the phase space. The mixing imphes and institutes average procedure. As a 
consequence, the information on separate phase trajectories is lost. This is equivalent to the 
irreversibility. The shortcoming of such an explanation is that the nature of averaging in a 
phase space because the dynamic equations do not contain the mechanism of coarse graining. 

Though the mechanism offered here is also based on mixing property, using a force, Fp, as 
evolution parameter, allows getting rid of explicit use of "coarse graining" idea. In contrast to 
directly "coarse graining" of the phase space, the transition to integration on impact parameters 
does not deform the nature of aspiration of the system to an equilibrium state. Here, the role 
of replacement of summation on integration on impact parameters consist in transition from 
discrete functions to continuous functions convenient for differential calculus. 
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